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ABSTRACT 



The Hubble Space Telescope has been used to obtain WFPC2 (V, /) photometry for 
a large sample of stars in the outer halo of the giant elliptical NGC 5128, at a distance 
4 Mpc beyond the Local Group. The globular cluster N5128-C44, at the center of 
the Planetary Camera field, is well enough resolved to permit the construction of a 
color-magnitude diagram (CMD) for it which covers the brightest two magnitudes 
of the giant branch. The CMD is consistent with that of a normal old, moderately 
low-metallicity ([Fe/H] ~ —1.3) globular cluster, distinctly more metal-poor than most 
of the field halo stars at the same projected location (which average [Fe/H] —0.5). 
This is the most distant globular cluster in which direct color-magnitude photometry 
has been achieved to date, and the first one belonging to a giant E galaxy. 



Subject headings: galaxies: elliptical and lenticular — galaxies: star clusters 



1. Introduction 

Globular clusters are typically found in much larger numbers within giant elliptical galaxies 
than in large spirals such as the Milky Way or M31, and in extreme cases such as cD galaxies the 
total population can exceed 20,000 clusters per galaxy {e.g. W. Harris 1991| ). It is conventionally 



assumed that the clusters in these remote galaxies must individually resemble the familiar Milky 
Way old-halo globulars in their ages, metallicities, and stellar compositions. But making direct 
tests of this fundamental assumption in different galaxies is a difficult task. For even the closest 
large galaxy, M31, it has become possible only recently to study the stellar compositions of its 
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old-halo clusters through direct color- magnitude photometry {e.g., Fusi Pecci et al. 199(^ ; Rich e\ 
al. 1996| ). The situation is worse for giant E galaxies, which are all far beyond the Local Group 



and in which individual clusters are seen as no more than faint, unresolved semi-stellar objects 
even under the best imaging conditions. The only reasonably direct tests of their nature are ones 
based on integrated light - color indices, spectral properties, line strengths - which verify that 
their characteristics at least roughly match those of the Milky Way globulars (see, e.g., pohen 
et al. 1997 for a comprehensive recent spectroscopic study of clusters in M87; or [Jablonka et al. 



1996 for spectral characteristics of clusters in NGC 5128). However, the information derived from 



integrated light must, by definition, be heavily smoothed over the composite stellar population. A 
color-magnitude diagram (CMD) for any such object, even with limited depth, would provide a 
uniquely direct test of their nature. 



2. Observations 



By far the closest giant elliptical galaxy is NGC 5128, the dominant member of the Centaurus 
group at d ~ 4 Mpc (the next nearest gE is NGC 3379 in the Leo group, which at d ~ 10 Mpc 
is a full two magnitudes more distant). Although its innermost few kiloparsecs are affected by 
the well known dust lanes and gas, the main body of the galaxy and its extensive halo are by all 



indications those of a basically normal giant elliptical {e.g., Graham 1979; Ebneter & Balick 1983 
Hui et al. 1995; Soria et al. 1996| ; Storchi-Bergmann et al. 1997). Since interactions and mergers 
of the type now happening in NGC 5128 are now recognized to be rather common occurrences 
in the ongoing histories of large galaxies, it is entirely reasonable to expect that its underlying 
stellar composition will be typical of gE's. In addition, its population of ~ 1500 globular clusters 



( G. Harris et al. 1984 ) has color and luminosity distributions that clearly resemble the clusters 
found in other, more remote gE's such as in Virgo and Fornax (see p. Harris et al. 1992| [hereafter 



HGHH], [Hesser et al. 1984| , and |H. Harris et al. 19881) . NGC 5128 unquestionably offers our best 
chance to investigate globular clusters in a giant E galaxy in the best possible depth and detail. 

In HST/WFPC2 program GO-5905, executed in two visits on 1997 Aug 16 and 24, we 
obtained deep images of a field in the outer halo of NGC 5128. The PCI CCD was centered on 
the globular cluster N5128-C44, which is located at a projected distance of Rgc = 18'32 from the 
center of NGC 5128, equivalent to 21.3(d/4) kpc. This location in the halo is far outside the 
core regions contaminated by the dust lanes and thus should be representative of the old-halo 
population of the elliptical (see HGHH for complete wide-field finder charts and a compilation of 
the relevant cluster data). C44 itself is a luminous, moderately low-metallicity cluster according 
to the photometry published in HGHH: at an integrated magnitude V = 18.48 {My = —9.87 for 
a distance modulus (m — M)y = 28.35; see below) it is similar to the luminous halo globulars in 
the Milky Way such as M3 or NGC 2419, while its (C — Ti) color yields an estimated metallicity 
[Fe/H] = —1.4 lb 0.25, quite similar to the mean metallicity of the Milky Way halo clusters and 
field stars. Lastly, the cluster radial velocity of 520 km s~^ is almost identical with the systemic 
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mean of ~ 540 km s ^ for its host galaxy (HGHH; |H. Harris et al. 198S ; Hesser et al. 198^ ). 



Om' WFPC2 observations consisted of 10 images totalling 12800 seconds in each of the two 
standard filters F606W {V) and F814W (/), with individual exposures sub-pixel shifted over 
five separate positions. The images were reregistered and combined to form composite deep 
(y, I) images free of cosmic rays and bad pixels. By simultaneously imaging both a metal-poor 
globular cluster and the halo field stars, our original intention was to gain the most direct possible 
comparison between their two mean metallicities, as well as to give us an additional consistency 
check on the metallicity scale from our {V — I) colors. Our photometry from the WF2,3,4 
images has been used to generate a color-magnitude diagram (CMD) for the galaxy halo, covering 
about 2.5 magnitudes of the old red giant branch (RGB) and containing more than 10,000 stars. 
That material is described in a companion paper along with a more detailed outline of the data 
reduction ( p. Harris et al. 1998 ); here, we present the results for the PCI field and the globular 
cluster. 

In the following discussion, we adopt a distance to NGC 5128 of 4.0 Mpc, based on planetary 
nebula luminosity functions (Hui et al. 1993), surface brightness fluctuations (Tonry & Schechter 



1990t Pui et al. 1993| ), and the luminosity of the RGB tip ( |G. Harris et al. 1998| ; ^oria et al. T996| ) 
all normalized to a contemporary Local Group distance scale in which ci(M31) = 770 kpc {e.g.^ 



van den Bergh 1995| ; [Fernley et al. 1998 ; and references cited there). The foreground reddening is 
adopted as E{B -V) = 0.11 and E{V - I) = 0.14, giving (m - M)y = 28.35. 



3. Data Reduction 

The image of the PC field around C44 is shown in Figure 1. The cluster subtends only a small 
angle (at d = 4 Mpc, typical globular cluster core and tidal radii of Tc ~ 2 pc, rt ~ 50 pc would 
have angular sizes of Od and 2'.'5, superimposed on the PCI image scale of 0'.'046 per pixel). Thus 
almost all of its stars are severely crowded and stellar photometry is restricted to the brighter 
stars on its outskirts. But the task is not hopeless. An excellent analogy of what we can expect 
to accomplish is to compare with the globular clusters in M31: in terms of spatial resolution and 
scale, imaging M31 from the ground under good seeing conditions of ~ 0'.'5 is closely comparable 
with imaging NGC 5128 (5 times more distant) with the O'.'l— resolution of HST. Successful 
ground-based studies of several M31 clusters were indeed achieved in the pre-HST era ( Pouture~et 



al. 1995; pieasley et al. 19*88 ; Christian fc Heasley 199l| ). In any one of these clusters, typically two 



or three dozen of the RGB stars could be isolated and measured, with representative photometric 
scatter of (j[V — I) = ±0.2 — 0.3 driven by the crowding and background nonuniformity. Although 
the precision of these studies was unavoidably low, it is important to note that the subsequent 
much deeper and more precise photometry carried out later with HST (see, e.g., Fusi Pecci et al\ 



1996; Rich et al. 1996|) verified that this pioneering ground-based photometry was able to produce 



systematically accurate colors and mean RGB loci. 
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We first extracted a 190 x 190 pixel subsection (8'.'74 x 8'.'74) from the PCI image, centered 
on C44 (Figure |l|). To eliminate as much of the unresolved light of the cluster as possible, we 
used the ellipse-fitting code within STSDAS to derive an isophotal contour model for the cluster, 
and subtracted this from the raw image. The result (second panel of Figure |l|) was then used for 
stellar photometry through the DAOPHOT II code ( [Stetson 1992 ). Three iterations of the usual 



sequence FIND/PHOT/ALLSTAR were employed, along with a very small psf-fitting radius (0.8 
fwhm), to find and extract the individual stars near the cluster. As a check on this procedure, 
we also tried removing the unresolved cluster light by a first pass of finding and removing stars, 
smoothing the remaining light by median filtering, subtracting that from the original image, and 



repeating the process (see Couture et al. 1995). The overall results for the stellar photometry were 
quite similar; however, the first (ellipse-fitting) procedure allowed us to reach slightly further in 
toward cluster center and to define the CMD a little more narrowly. 

Zeropoints and color terms for the conversion of instrumental magnitudes into {V, I) followed 



the prescriptions of Holtzman et al. (1995) and are discussed completely in G. Harris et al. (1998). 



4. The Color-Magnitude Diagram 

The resulting CMD is shown in Figure |2|, in which all measured data points are plotted 
without any selection criteria applied. The stars closest to the cluster (r < 1'.'4 ~ 30 px, which 
we will call the "inner zone") are systematically bluer than the surrounding field-halo stars, 
suggesting that this particular cluster has distinctly lower metallicity than the average of the halo 
stars. To define a cleaner CMD for the cluster, we culled the data in two additional ways: (a) 
stars with very poor goodness-of-fit values (x > 4, although the great majority of stars had x < 2) 
were eliminated; and (b) individual stars were removed from the inner-zone CMD in proportion 
to their presence in the outer (r > 30 px) zone. For this latter step, we divided the CMD into 
a grid of boxes of size Al ~ 0.5, A{V — I) ~ 0.2. Then, noting that the area of the outer zone 
(70.0 arcsec^) is ~ 11 times larger than that of the inner zone (6.4 arcsec^), we removed one star 
from the inner-zone CMD in any box which also contained ~ 10 it 3 stars from the outer zone. 
(This star-by-star extraction will actually remove too few stars for 25.5, where the detection 
completeness is much lower in the highly crowded cluster than in the surrounding field. Thus we 
expect the fainter parts of the cluster CMD to remain partially contaminated, but by only two or 
three stars at worst.) 

The "cleaned" CMD for the cluster is shown in Figure ^. The 45 remaining points now 
scatter more closely around the fiducial line for [Fe/H] = —1.3 over the 2.5-magnitude run of the 
data, particularly at the bright RGB tip where the stellar photometry should be the most reliable. 
This metallicity estimate is uncertain by perhaps ~ ±0.3 dex taking into consideration both the 
random photometric scatter and the uncertainties in reddening and distance moduli; but it is 
entirely consistent with the previous metallicity calibration [Fe/H]c_r^ = — 1.4 from its integrated 
colors (HGHH). The typical color uncertainty a{V — /) ~ ±0.15 — 0.2 per star, along with some 
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residual contamination by field stars, is enough to be responsible for the remaining scatter we see 
in the CMD. 



5. Radial Profile and Structural Parameters 

Cluster C44 is sufficiently extended on the PCI image that we can make useful estimates of 
its structural parameters. The smooth profile of the cluster light generated from our stsdas. ellipse 
and . bmodel fits is shown in Figure Q . Both the V and / surface intensity plots indicate that 
the core of the cluster is significantly more extended than the PSF profile, and at large radii, the 
light can be traced out past r ~ 2" (40 px). The profile shape is also nearly circular and fairly 
uniform with radius, as shown in Figure |5|; the eccentricity of the model fit goes through modest 
fluctuations about its mean value (e) ~ 0.06. 



Comparison with representative King (1966) model profile curves, also shown in Figure ^ 



indicates that the cluster is adequately matched by models with central potential parameters in 
the range Wq = 7 =b 1. A very rough estimate of the core radius of the cluster from the model 
fits is Tc ~ O'.'OS ~ 1.5 pc, after subtraction in quadrature of the PSF "core" rc^psF — 0'.'04. The 
half-mass radius, which is larger and thus more well determined, is at r^ ~ 0'.'27 ~ 5 pc, and 
the central concentration index of the fitted curve is c = log{rt/rc) = 1.6 it 0.3. Its estimated 
tidal radius is then rt ~ 3" ~ 60 pc. These values all lie comfortably within the typical range for 
globular clusters in the outer halo of the Milky Way, which have more extended cores and larger 
tidal radii than those in the inner bulge {e.g., [Trager et al. 1995 ; van den Bergh et al. 1991). 



In summary, this study clearly favors the important "null hypothesis" that N5128-C44 is 
an old, metal-poor globular cluster closely resembling the familiar ones in the Milky Way in 
every way we can test it: metallicity, integrated color, color-magnitude diagram, and King-like 
structural parameters. This is the first time that direct stellar photometry has been used to test 
this hypothesis in any globular cluster in a giant E galaxy. Our results support the contention that 
we are genuinely dealing with the same fundamental type of object when we perform comparisons 
between old-halo globular cluster populations in gE galaxies and those in the Local Group. We 
suggest that the photometric analysis of globular clusters in other galaxies well beyond the Local 
Group is within reach, and will become routine with the deployment of more advanced imaging 
tools. 
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Fig. 1.— (a) First panel: The central 190 x 190 pixels (8'.'7 x 8'.'7) of the Planetary Camera WFPC2 
field, centered on globular cluster C44 in the halo of NGC 5128. The image is a 12800-second 
sum of 10 exposures in V (F606W). The smooth unresolved light of the cluster has been partially 
subtracted to emphasize the resolution into stars of its outer envelope. The inscribed circle denotes 
the radius (30 px = 1'.'4) of the inner zone used to define the cluster CMD. (b) Second panel: The 
same PCI field after removal of elliptical isophotes. 
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Fig. 2. — Color-magnitude data for all measured stars in the central PC field. Solid dots (left 
panel) denote stars within the inner 30-px (l'.'4) circle; crosses (right panel) denote objects outside 
that circle, which are mostly RGB stars in the halo of NGC 5128. The four lines superimposed 
on the diagram are fiducial lines for Milky Way globular clusters of four different metallicities (Da 
Costa fc Armandroff 1990| ; puarnieri et al. 1998| ); from left to right, they are M15 ([Fe/H] = -2.2), 



NGC 1851 (-1.3), 47 Tuc (-0.7), and NGC 6553 (-0.25). The lines have been placed assuming 
an intrinsic distance modulus for NGC 5128 of (m — M)o = 28.0 and reddening E{V — I) = 0.14 
(see text). 
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Fig. 3. — Statistically cleaned CMD for the globular cluster C44. Stars with poor PSF fits (x > 4) 
have been rejected, and background contamination in proportion to the area of the inner zone has 
been removed as described in the text. The fiducial cluster sequences for four different metallicities 
are the same as in Fig. 2. Photometric measurement uncertainties from the ALLSTAR fitting 
algorithms are indicated for three different magnitude levels by the error bars at right. C44 appears 
to have a moderately low metallicity [Fe/H] ~ —1.3 based on the color of its RGB, and in close 
agreement with the metallicity from its integrated (C — Ti) color index. 
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Fig. 4. — Radial profile of the cluster light, plotted as surface brightness (magnitudes per unit area 
relative to the central intensity) against log A (semi-major axis) in arcseconds. Solid dots denote 
the surface intensity in V of the smoothed isophotal model of the cluster; open circles denote the 
intensity profile in /. The dashed line at left shows the profile of the stellar point-spread function. 
Solid lines superimposed on the data points are profiles for four King models with different values 
of the central potential parameter Wq (curves of increasing central concentration go from left to 
right). We estimate the cluster to have a core radius Vc — C.'OS (about twice as large as the PSF 
core) and tidal radius rt ^ 3". 
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Fig. 5. — Eccentricity e of the cluster profile, as a function of semi-major axis A. The e— values are 
deduced from the elliptical isophote fit to the cluster light described in the text. The cluster has a 
small and nearly uniform eccentricity (e) 2± 0.06. 



